f{R) cosmology with torsion 
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/(-R)-gravity with geometric torsion (not related to any spin fluid) is considered in a cosmological context. 
We derive the field equations in vacuum and in presence of perfect-fluid matter and discuss the related cosmo- 
logical models. Torsion vanishes in vacuum for almost all arbitrary functions f{R) leading to standard General 
Relativity. Only for f(R) — R^, torsion gives contribution in the vacuum leading to an accelerated behavior . 
When material sources are considered, we find that the torsion tensor is different from zero even with spinless 
material sources. This tensor is related to the logarithmic derivative of f'{R), which can be expressed also as a 
nonlinear function of the trace of the matter energy-momentum tensor E^i, . We show that the resulting equa- 
tions for the metric can always be arranged to yield effective Einstein equations. When the homogeneous and 
isotropic cosmological models are considered, terms originated by torsion can lead to accelerated expansion. 
This means that, in f(R) gravity, torsion can be a geometric source for acceleration. 
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I. INTRODUCTION 



O ' ACDM has recently assumed the role of a new Cosmological Standard Model giving a coherent picture of the today observed 
universe 01]. Although being the best fit to a wide range of data, it suffers of several theoretical shortcomings |2] so it fails 
in tracking cosmic dynamics at every redshift and fails in according observational cosmology to some fundamental theory of 
physical interactions. Among the defects of this model, there is the lack of final probes, at fundamental level, for dark energy and 
dark matter candidates (which should be the 95% of the energy-matter content of the universe!) which frustrates the possibility 
to reduce ACDM to some self-consistent scheme, despite of the fact that it is a fair "snapshot" of the present status of the 
J> universe. These facts motivate the search for other models, among which alternative theories of gravity that should reproduce 
0^ the successes of ACDM but should be more appropriate in describing the cosmological dynamics 121 01 • 

In particular, the large part of dark energy models relies on the implicit assumption that Einstein's General Relativity (GR) 
is the correct theory of gravity indeed. Nevertheless, its validity on large astrophysical and cosmological scales has never been 
tested but only assumed Jst], and it is therefore conceivable that both cosmic speed up and missing matter, respectively the dark 
energy and the dark matter, are nothing else but signals of a breakdown of GR at large scales. In other words, GR could fail 
in giving self-consistent pictures both at ultraviolet scales (early universe) and at infrared scales (late universe) also it is fairly 
working at Solar System scales and in the weak field regime. 

Staring from these considerations, a different possibility could be to better investigate the gravitational sector and consider, as 
source of the field equations, only the observed (and probed at fundamental level) baryonic matter and radiation (photons and 
neutrinos). A choice could be to take into account generic functions f{R) of the Ricci scalar R. The goal should be to match 
observational data without considering exotic dark ingredients, unless these are found bymeans of experiments at fundamental 
level This is the underlying philosophy of the so-called f{R)-gravity (see e.g. iflMI^)- 

These theories are receiving much attention in cosmology, since they are able to give rise to the accelerating expansion ff^ 
and it is possible to demonstrate that they play a major role also at astrophysical scales (for recent comprehensive reviews see 
[Toll ). In fact, modifying the gravity Lagrangian affects the gravitational potential in the low energy limit. Provided that 
the modified potential reduces to the Newtonian one on the Solar System scale, this implication could represent an intriguing 
opportunity rather than a shortcoming for f{R) theories. In fact, a corrected gravitational potential could offer the possibility to 
fit galaxy rotation curves without the need of huge amounts of dark matter llllllT2llT3lll4l1T5iil6lfl7ll . In addition, it is possible 
to work out a formal analogy between the corrections to the Newtonian potential and the usually adopted galaxy halo models 
which allow to reproduce dynamics and observations without dark matter 1 13]. 

However, extending the gravitational Lagrangian could give rise to several problems. These theories could have instabilities 
ifTsll . ghost -like behaviors [19], and they have to be matched with the low energy limit experiments which quite fairly test GR. 

In summary, it seems that the paradigm to adopt /(i?)-gravity leads to interesting results at cosmological, galactic and Solar 
System scales but, up to now, no definite physical criterion has been found to select the final f{R) theory (or a class of theories) 
capable of matching the data at all scales. Interesting results have been achieved in this line of thinking [2CI, .21j 22, 23, 24] but 
the approaches are all phenomenological and are not based on some fundamental principle as the conservation or the invariance 
of some quantity or some intrinsic symmetry of the theory. In any case, some results are emerging in this direction and the 
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pres ence of Noether symmetries into dynamics seems to play an important role in selecting physically interesting f{R) theories 

imm. 

Besides selection criteria, the "full" geometric sector of f{R) gravity has to be investigated considering also the role of 
torsion. Such an "ingredient" has been firstly considered by Cartan and then by Sciama and Kibble in order to deal with spin in 
GR (see IztIi for a review). Being spin as fundamental as the mass of the particles, torsion was introduced in order to complete 
the scheme that mass(energy) are the source of curvature and spin is the source of torsion. Unfortunately, torsion in the context 
of GR, seems not to produce models with observable effects since the gravitational coupling is extremely weak in all the torsion 
phenomena and only in the very early universe its effect could have been significant. 

However, it has been proven that spin is not the only source for torsion. As a matter of facts, torsion can be decomposed in 
three irreducible tensors, with different properties. In |28], a systematic classification of these different types of torsion and their 
possible sources was discussed. 

In two recent papers, the role of torsion in f{R) gravity has been considered in the framework of metric-affine formalism 
ll29i I30I1 . The field equations have been discussed in empty space and in presence of perfect fluid matter taking into account 
the analogy with the Palatini formalism llsill . As a result, the extra curvature and torsion degrees of freedom can be dealt as an 
effective scalar field of fully geometric origin {2^. 

From a more formal viewpoint, f{R) gravity with torsion can be studied in the framework of the J^-bundles formalism llsoll . 
Such an approach is particularly useful since the components of the torsion and curvature tensors can be chosen as fiber 
coordinates on the bundles and then the symmetries and the conservation laws of the theory can be easily achieved. Also in 
this case, field equations of /(i?)-gravity have been studied in empty space and in presence of various forms of matter as Dirac 
fields, Yang-Mills fields and spin perfect fluid. Such fields enlarge the jet bundles framework and characterize the dynamics. 

In this paper, we discuss the cosmological applications of /(_R)-gravity with torsion, considering the possibility that the whole 
dark side of the universe (dark matter and dark energy) could be geometrically interpreted by curvature and torsion. Exarnples in 
which repulsive gravity and clustered structures could be implemented considering torsion are present in literature ||27[ l32i [33ll 
but, in that cases, GR has been adopted and the whole dark sector has not been addressed. 

The layout of the paper is the following. In Sec. II, we give the definitions and conventions for torsion tensor and construct 
the in the main geometrical quantities in spacetime. Field equations in vacuum are derived in Sec. III. In Sec. IV, we derive 
the cosmological equations in vacuum and find out some interesting solutions where cosmological term is given by the trace 
of torsion tensor Field equations, in presence of standard fluid matter, are discussed in Sec.V, while, in Sec. VI, also spin is 
considered. A relevant discussion is devoted to the fact that torsion contribution can be dealt under the standard of a scalar field. 
In Sec. VII, we consider weak and strong energy conditions in f{R) gravity with torsion. Cosmological models in f{R) gravity 
with torsion and exact solutions are discussed in Sec. VIII. Discussion and conclusions are drawn in Sec. IX. 



II. TORSION TENSOR AND INVARIANT QUANTITIES IN U4, 



A 4D-differential manifold equipped with torsion is defined in a U4, space, while standard torsionless Riemannian manifolds 
are defined in V4. In literature, there are several definitions for torsion and quantities related to it. For a summary, see ll28ll . The 
conventions adopted in this paper are those in [|27| . 

The torsion tensor 5^^'^ is defined as the antisymmetric part of connection in a coordinate basis 

a connection with torsion and metric-compatible (Va5^i/ — 0) has the form 

^ fiiy {m'^ } ^ ^fiiy ^i' fj, ^ ^ fj,u {^i^ } ^ (^) 

where {^^} are the Christoffel symbols and 

liu ' U /J ^ flu V"^/ 

is the contorsion tensor. Another combination of torsion tensor, often used in the calculations, is the modified torsion defined by 
the following relation 

= S^jj + &^Si, — Sj^Sfi (4) 

where 



(5) 
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is the torsion trace-vector. 

Here, we shall consider a vectorial torsion |28] of the form 

S^/ = A[,S^] . (6) 

Its trace is 

3 

so the vectorial torsion takes the form 

In the following, we will find convenient to express the field equations in terms of one of these quantities. To this purpose, let us 
give here the relations among torsion, modified torsion, contorsion, and their respective traces. From Eq.©, it is 



and then 



T^^-2S^ or S^ = -^T^ (10) 



where the modified torsion can be expressed in the following form 



t^.^ = It[X]- (11) 

The Riemann tensor can be decomposed in a part depending only on the Christoffel symbols and a part which contains the 
contorsion tensor. Using Eq. lfTTT i in Eq.©, we obtain 



K^.^ = l{T^.^-T,\+T\,) . (12) 



Finally, substituting (fTTT i in (fT2] i. we get 



K^,^ = i - S^^T,) (13) 



and 

= g^'^K^/ = T^. (14) 
Furthermore, starting from the standard definition of the Riemann tensor 

one can insert the decomposition of the connection in the Christoffel and contorsion parts 

Kf3 = {^,}-Kp"- (16) 

In this way, the Riemann tensor results decomposed in a first term obtained by the Christoffel connection and its derivative and 
a second term given by the contorsion and its derivative. 

= d, ({ - i^./) - 9. - ^,/) 

+ {{^0} - ^./) - ^M.") - ({,^} - ^,/) - Ka") ■ (17) 

Using the covariant derivative without contorsion V, the previous expression takes the form 
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The corresponding Ricci tensor and the curvature scalar are respectively 



and 



(19) 



R ^ g^^^Rp, = Rp,{{}) + g^'^'V.X,/ - g^N^K^^" + g'^^K^^/K^," - g^^^K^/K^," . 
If torsion is in the vectorial form, the Ricci tensor and the scalar curvature take the forms 



(20) 



(21) 



R = i?({}) - 2VaT^ - -T^Tx , 



which wiU be widely used in the following discussion. 



(22) 



III. FIELD EQUATIONS IN THE VACUUM 
Let us start our considerations deriving the field equations in vacuum. The action for the gravitational part is 

^ = y V^gf{R)d^x. 



(23) 



It is instructive to report the variational principle step by step in order to put in evidence the differences with respect to a variation 
in V4. It is 



5A 



I 



f'{R)R^. - -g^..f{R) 



5g 



^"'d''x + j ^f\R)g^''5R^,d^x 



(24) 



(25) 



/ 



f'{R)R^, - -g^.f{R) 



5g^^^d''x + j sf^gnR)g^'' [2V[.5r-], 



+2S,^>^5T^,^]d^x. 
In a space equipped with torsion, the following property holds 

d{P5Q) = PVSQ + {VP)SQ , 



(26) 



(27) 



where PSQ is a vectorial density and = + 25^^"^. If the variation SQ is zero on the boundary or at infinity, we find 

P\/6Q = - {VP)dQ . (28) 



The variation assumes the form 



6A 



nR)R^. - -g^.f{R) 



.d^x 



(29) 
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Giving the explicit expression for V, we obtain 

8A^ 



(30) 



where T^^ is the above modified torsion tensor. Finally, the variation of the connection, expressed in terms of the metric and 
the contorsion, is 



1 



By substituting Eq.dSTI) in ( l30b . we obtain 



(31) 



(32) 



In general, f{R) theories, in metric formalism, present fourth order terms in the field equations. In this case, such terms are 
absorbed in the torsion components and then 



The divergence in Eq.(l33Tl is zero because of Eq.(l32b and then the field equations in vacuum are 
and 



(33) 



(34) 



f'iR}T./+ 2 dxif'iR)) {S'.Sii - S^6P,) = . 



Taking into account the trace of Eq.(l34ll 
and substituting it in Eq. (|34] |. it follows 



/'(i?)i?-2/(i?) = 0, 



/'(i?) ( R^. - -g^,R ] = 



(35) 



(36) 



(37) 



The set of Eqs.(l36]l and ( |37] ) is equivalent to Eq.(l34b. 

In general, Eq.([36]l is an algebraic or transcendent equation for R which is satisfied if and only if i? is a constant. This means, 
from Eq.(l35]), that torsion is identically zero and then we have. The resulting spacetime may be non-trivial, say a Schwarzschild- 
de Sitter. The trivial case may be de Sitter or anti-de Sitter, according to the sign of A, i.e. 



R^v ^9f2iy ■ 



(38) 



Only for f{R) = aR^, Eq.(l36]l is an identity and R is not necessarily a constant. In this case, torsion is different from zero in 
the vacuum. By substituting f{R) = aR? in Eq.(l35ll. we have 



which can be reduced to 



" 2 R 



(39) 



(40) 
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This equation, together with 

R^. - \9f..R = (41) 

is a system of equations for the gravitational field. 

We can decompose the curvature terms in Eqs. (|4TI) and ( |40l ) in a part given by the Christoffel symbols and a part depending 
on torsion. The resulting equation is 

RM) - ^.9p.i?({}) - p.T^ + ^5p.V^T, + ^T^n - ^9,uT^Tx - (42) 

where the symbols V and i?({}) represent, respectively, the covariant derivative and the scalar curvature given by the Christoffel 
symbols. Similarly, the Ricci scalar R results 

R = R[{}) - 2VaT^ - ^T^Ta . (43) 
On the other hand, the vector component of torsion satisfies the following differential equation 

^ a, {r{{}) - 2^xT^ - ^^aT^^ = (^R{{}) - 2VAr^ - ^TaT^^ (44) 

which is a propagation equation for T^.. Finally, applying the contracted Bianchi identities, adding and subtracting _R({})/4 in 
Eq. (l42] i and taking the covariant divergence with respect to the V derivative, it follows that 

^^''^V.VaT^ - l^AT^T^ - \d^R - ^9^V,T^ + ^9^(TaT^) = , (45) 
which we will be useful for the discussion below. 



IV. COSMOLOGY WITH TORSION IN VACUUM 



In order to develop our cosmological considerations, let us take into account a Friedmann-Robertson-Walker metric of the 
type 



The scalar (Christoffel) curvature in such a metric is 

i?({}) = -6 



[dx^ + dy^ + dz^] . 



(46) 



d 




k ' 






- + 






= -6 (^iJ + 2H^ H 




a 






a2 J 



(47) 



Being the system dependent only on t, from Eq.(l44li. only the Tq = T component of is different from zero. Eqs.(l42]l and 
( |47] i, after some algebraic manipulations, give the following cosmological equations 



~ {%H + 12i?2 + 64r + 2r + &HT + ^T^ 
2 at \ 3 



+ 12H^ + + 2t + &HT +rT'^ ]T , (48) 



'da k 



-t +-HT +-T'^ . 
3 3 9 



a a a 

which, in terms of the Hubble parameter H, can be reduced respectively to 

k k 
3H + \2HH - 6^H + T + HT + 3HT = AH'^T + 2^T + 2HT^ 



rp6 



(49) 



(50) 



and 
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H -^ = -lf + \HT+lT^ (51) 
3 3 9 



The Bianchi identities ( |45T l gives 



.... . . .. . h 

3T - 2TT + \2HT - UH^T - 3HT - 3HT - AHT^ + 9H + 36HH - 9H^ = (52) 

while the spatial components reduce to an algebraic identity. Let us find a solution for the system ( |50b , ( |5T] i ( fSTT i in the physically 
interesting case fc corresponding to a spatially flat universe. Substituting H in ( l52b , it reduces to 

^^^^H+1;tX ^1(h+It]\, (53) 

and then we have two cases 



dt \ 3 J 3 V 3 



H+\t = 0, (54) 



'h+]-T] = \ (h+]-T]T, (55) 



In the first case, we have the solution 



In the second case, we obtain 



3 

, , -T , , , 

dt\ 3 / 3 V 3 



a(<) = aoexp( --Tt ) (56) 



F = i/oexp( It") -| (57) 



which gives a cosmological expansion driven by torsion; 

ait) = gq exp 



T , /T 
--^i + ^o exp l-t 



(58) 



being _ffo> ™d arbitrary integration constants. This situation is particularly interesting in view of considering vectorial 
torsion as the geometrical source of cosmological accelerated expansion also without assuming an additional spin fluid Issll . It 
is important to point out at this point that solutions ( |56] ) and ( l58l ) are derived in the particular case of T =constant but they can 
be derived also for T = T{t). In this general case, T has to be substituted with T = J T{t)dt. These results are clear indications 
that geometric torsion can be an effective source for inflation and, in general, accelerated expansion. 

V. FIELD EQUATIONS IN PRESENCE OF PERFECT-FLUID MATTER 

A more realistic situation is considering, in the action ( |23] |. the presence of a perfect fluid matter Lagrangian density. We have 

A= I V^ifiR) + Cm]d^x . (59) 



By the sake of simplicity, let us consider first £,„ not containing torsion terms. The corresponding matter fluid energy-momentum 
tensor is E^'' and E — g^^Yj^^'^ its trace. The field equations are then 



- -9^,,f{R) = S^, (60) 



and 



nR)V + \ dxwm {s'.s'a - s^sf,) = . (61) 
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The trace of Eq.(l60]l 

/'(i?)i?-2/(i?) = S, (62) 
gives a non-linear relation between the curvature scalar R and S which is 

R = ^(S) , (63) 

and then Eq.dSOb and ( |6TI ) can be rewritten as 

1 1 / 1 \ 



4^/^"^ f'{F{^)) 

V--l'-J^^m-^y.)- (65) 

It is worth noting that the equation for torsion is now an "algebraic" expression. Torsion is present as a sort of gravitational 
coupling being Eq.(|64b can be written in an Einstein form by adding and subtracting —\g^j,vR and then using 

Eq.(l63l). We obtain 

1 If 1 \ 1 , , 

Ri.u - i^9i.uR = f'{F{Y.)) I ^^"^ ~ i^^''^ ) ~ i^f""^^^' ■ ^^^^ 

By using Eqs.(l42]i and ( l43T l. we can decompose the Ricci tensor and the curvature scalar in their Christoffel and the torsion 
dependent terms 

RM) - l9^..R{{}) - f'(^F{j:)) i^^" " l^'^"^) " + ^^^^'^ 



By taking the trace of Eq. 



where we define the auxiliary scalar field 



we obtain 



-gffM-VAT^ - -T^T, - -g^uTxT\ (67) 



2 ip 



V = nF{E)) , (69) 



1 31 31 1/S \ 

ip 2 (p'' 4(^^ 4 \(p J 

+ -Vf,d^p- -gf_cuOtp. (70) 

The RHS of this equation can be rearranged as the sum of the energy-momentum tensors of standard perfect fluid matter and 
of a scalar field being 



1e -J- 

(p 2p^ 



2 ~ 2 ~ ' 

-pVfj.dup}^' -g^^pUp 



(71) 



where the scalar field effective potential is given by 



V{p) = i {p'{f')-\v) + V'i^-^d/')"'!^))) (72) 

and the operator □ = Va9'^ is defined, as above, by the covariant derivative without contorsion. This result shows that our 
approach is fully equivalent to an effective theory with perfect fluid matter and a scalar field with a geometrical interpretation 
given by Eq.(|69]l and dynamics given by the self-interaction potential ( |72] |. We obtained this picture without adopting conformal 
transformation but only splitting torsion from Christoffel contributions in the connection. In other words, the result strictly de- 
pends on the intrinsic non-linearities of the gravitation interacting with perfect-fluid matter assuming a gravitational Lagrangian 
density of the form \J—gf{R) defined on a Ui manifold. 
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VI. FIELD EQUATIONS IN PRESENCE OF MATTER AND SPIN 

As it is well known, spin can be one of the sources of torsion. In 1923, Cartan showed that intrinsic angular momentum 
(after called the spin) could have an important role in a geometric theory of space-time like General Relativity ll34ll . He showed 
that spin could originate torsion. This idea was considered also by Sciama [35] and Kibble |36] and, in more recently, by Hehl 
and his coworkers |27] and Trautmann and his coworkers |37]. The path followed in this stream of research was to extend GR 
to a theory (Einstein-Cartan-Sciama-Kibble, ECKS, theory) in which spin is the source of torsion. In the limit of zero spin 
distribution, the theory reduces to standard GR. 

In the ECSK theory, spin and torsion are related by an algebraic equation and torsion does not propagate. In this case, only 
considering a distribution of aligned spins, described by a spin-density tensor t^^ it is possible to deal with torsion as a fluid. It 
is then necessary not only to have fluid of particles with spin, but also the average spin alignment has to be different from zero, in 
order to produce torsion. To deal with spin-torsion cosmology in the ECKS approach, it is necessary to consider a cosmological 
fluid where spin alignment is produced, for example, by huge magnetic fields. Such field are difficult to be found in standard 
cosmological situations as discussed in several works in literature. 

In 1 28], it was discussed in details that spin is not the only source of torsion and that the torsion tensor can be decomposed 
in three irreducible tensors: i) the totally antisymmetric torsion, ii) the vectorial torsion and in) the traceless torsion. It was 
shown, for example, that a spin field generates a totally antisymmetric torsion, while a classical spin fluid, as that introduced by 
Mathisson |38] and Weyssenhoff |39], generates a traceless torsion. 

In the previous section of this paper, we have seen that the vectorial part of torsion can be generated by the logarithmic 
derivative of a given scalar field of geometric origin. 

Considering also spin in the context of this paper leads us to have a combined effect of different sources of torsion. In other 
words, one has to take into account the vectorial torsion, obtained by the covariant derivative of ip, and the classical spin average 
effect. If in the action ( |59l ), we take into account a material Lagrangian £,„ with spin, the field equations have to be modified in 
the following way 

SV^f{R) ^ _f'(^ji-^jit^i' + -gi^''f[R) + VxiU^""^ - U"^" + C/^^'') = i === ^ (73) 

Sg^iu 2 2 ^ Sg^^y 

U,/ EE ^^^^ = g'^dxifiR)) (S^iS^, - 5X) + V'{R)g^''V = ^J§^ ^ ■ ^^4) 

oK^p yjg OJ^p 

Such equations result more complicated than Eqs.(l60t. ( l6Tl l. In this paper, we will not consider anymore the presence of the spin 
fluid as one of the source of torsion. The remarks of this section have been necessary in order to point out that torsion can have 
also a purely geometric origin. In this case, despite of the spin fluid which gives only an algebraic relation, the torsion dynamical 
field related to the geometric part can greatly contribute to the cosmic dynamics being a very natural source for cosmic speed up. 

VII. THE WEAK AND STRONG ENERGY CONDITIONS IN f[R) GRAVITY WITH TORSION 

As we said, the f{R) gravitational field equations with torsion can be expressed in the equivalent Einstein-like form 



G ni^ ^liU p. o^P-^ (75) 

Lp 2 



where a^j^^, is the effective scalar field energy-momentum tensor in Eq.dTTTi. Such a term has a purely geometric origin and takes 
in account the nonminimal coupling between standard fluid matter and /(i?) gravity. It is important to stress that it gives a 
negative contribution in RHS of Eq. ( |75] ) and it is possible to show that it can naturally give rise to the observed accelerated 
behavior (due to ~ 0.7) and the clustering properties (due to Vl„i ~ 0.3) of the standard ACDM. 

In order to consider the properties of a^i, and its relevance for the cosmic evolution, let us recast it as a (perfect) fluid energy- 
momentum tensor. Following this prescription, after considering Eq.(|68j, let us define 



= (76) 
which plays the role of a. four-velocity. We are assuming that T'^ is a timelike-vector Then the energy density is 
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(f 2 (f^ if 2 Lp-^ \2 3 



(77) 



and the pressure 



P 1 3 1 



(78) 



For our purposes, it is important to consider the weak and strong energy conditions [40]. The weak energy condition implies that 



P > 



Q> 



3 1 
2^ 



1 



(79) 



where g = Yj^yU^V^ and p — Y.^i,{g^'' + U^U'^) are respectively the energy density and the pressure for the standard-fluid 
matter source. 

The strong energy condition holds, if the RHS tensor of Eq. ( fTST i satisfies the condition 



3 


1 






3 




2 








2 





(p + 3P) > 



(80) 



corresponding to 



3 1 

e + 3p> -- 

2 ip 



dxpd^ip ~ V{p) - ^ {^p> Up> + (^(9A¥'9V)"'5^'^aVV^9,(p 



(81) 



When this condition is satisfied, the universe undergoes a decelerated cosmological expansion, viceversa an accelerated expan- 
sion is expected when it is violated. These inequalities strictly depend on the form of the function /(i?) and on the nonminimal 
coupling of f{R) with the standard matter given by S (see Eqs.(|62l) and (|69l)). 

In conclusion, we can say that introducing torsion of geometric origin in f{R) gravity can lead to an accelerated behavior 
of the universe due to a repulsive nonlinear interaction of the (baryon) matter with itself. A qualitative discussion of this point 

^ E V{w) 

is given in [29]. There is discussed also the possibility that comparing the effective Se// — — and Ve// = — — can lead to 

comparable values for 51m and JIa. This could be a straightforward way to solve the so-called coincidence problem d |^ [H 0] • 



VIII. THE COSMOLOGICAL SOLUTIONS IN PRESENCE OF PERFECT FLUID MATTER 

In order to determine the cosmological equations in presence of perfect-fluid matter, we can arrange Eq. dTOb in the form 



^ 4 V<^ ^ ' ) ^ 2 ip^ 



ip 2 LP ^ " ^'-^ ip 2 ^'^ ip 
We consider the components of the Ricci tensor 



5';-9'''n^— + A^5^^. (82) 



and 



Rl ^ iso - ±E + -^-"£^ + 11 + ^H^ (83) 

^ip°Aip 4,^'2ip^ 2ip 2ip ^ ' 



^ ip ^ 4,Lp A ^ ' 2ip 2 ip ^ ' 
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In the homogeneous and isotropic FRW universe, we have 

n 

(85) 

(86) 

(87) 
(88) 
(89) 

where < 7 < 1 is the adiabatic index of the equation of state. The cosmological equations are then 



i?0 = 


-3-, 






Rl = 


-il- 




2k\ 


SO = 


Q, 






= 


-P = - 


-(7- 




E = 


g-3p 


= (4- 


- 37)e, 



a 3 70 1,, i ilP' 3(3 3 (3 
3- = -^ + -F I] + + (90) 



7 e 



1 



4 4 

By composing these two equations to eliminate the a /a, one obtains 



2 



4 12 
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which, solved with respect to H, gives 
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In order to derive the evolution of the density g, let us consider the following equation, derived by combining Eqs.(|90ll and ( |9T1 ), 
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This equation simplifies to 
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A further simplification comes from writing, explicitly in Eq. 
equation 



the expressions for A^, and H. We finally find the evolution 
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These equations can be solved by quadrature. The scale factor a{t) is obtained by solving Eq. 
As an example, let us assume f{R) oc i?" and consider a perfect fluid with equation of state p = (7 — l)g in a universe with 
k = 0. After some calculations we find, from Eqs. ( |63] ) and ( |69] l. 
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and 
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with 
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Solving Eq.(|97l). we find 
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The explicit expression for (p(t) is immediately found and 



(C2C|(n-2)2) 1 



72 (t - tor ■ 

The evolution of a{t) is obtained by integrating Eq.(|93]l. We have solutions for n > 4, that is 
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and then 



a{t) = Ai(t ~ to)" where a = n — 1 + e 



27 {t - to) ' 

{CiC2{n-2)) 



27 



(101) 



(102) 



(103) 



(104) 



(105) 



This is a power law expansion depending on the parameters 7 and n. Accelerated/decelerated behaviors are easily achieved by 
discussing the values of these parameters with respect to the Bianchi identities. As an example, let us consider the case n = 4. 
The decelerated behavior is obtained for 



0<a = 3 + £^^ < 1 



(106) 



otherwise the universe is accelerating. The transition from the decelerated to the accelerated phases is achieved for a ~ 1 and 
then the redshift has to evolve as 2: t~^. Following 1 21], it is easy to see that the condition [1 + r^e//] — 0-2 to recover ACDM 
prescriptions is easily recovered depending on the couple of parameters {n, 7}. 



IX. DISCUSSION AND CONCLUSIONS 



Several issues from cosmology and astrophysics are telling us that the Einstein General Relativity should be revised in order 
to avoid shortcomings as dark energy and dark matter which, up to now, have not been probed at any fundamental level, but 
manifest their presence at large scales. Despite of the lack of final experimental evidences, the dark side of the universe should 
constitute almost the 95 % of the whole matter-energy content. 

This situation is extremely disturbing and a more "economic" way to solve the problem could be to revise the geometric part 
of the gravitational interaction extending the GR. In f{R) gravity, this point of view is considered and, in general, the good 
results of the Einstein theory are preserved trying, at the same time, to encompass the observational data in some self-consistent 
scheme |@]. 

In this paper, we have pushed forward this approach taking also into account the role of torsion in f{R) cosmology. With 
respect to the ECKS theory, where torsion couples algebraically to the spin and matter without dynamics, we have put in evidence 
the fact that torsion fields of purely geometric origin can perfectly mimic the role of scalar fields in cosmological evolution. Such 
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a feature naturally gives rise to accelerated behaviors as in dark energy models. In general, torsion has no effect in vacuum but, 
for R^, it is a source for the field equations leading to de Sitter-like expansions. 

In presence of standard fluid matter, torsion and curvature degrees of freedom give rise to nonminimal couplings which deter- 
mine cosmological evolution. This depends on the form of f{R). Also in these cases, accelerated behaviors are easily achieved. 
However, also in presence of torsion, /(i?)-gravity models must satisfy some viability conditions in order to describe cosmic 
acceleration. Such conditions are summarized [? ]. The cosmological models proposed here fully satisfy these compatibility 
conditions since the presence of torsion acts as a further additive fluid with the only effect to move the extremum of the effective 
potential which, following the notation in [? ], lies at the GR value R ^ [p — 2>p). 

As a final comment, it is worth noticing that we have not used any conformal transformation so the issue to choose between 
the Einstein and the Jordan frame is avoided. Besides, torsion naturally results as a fundamental scalar field whose origin is 
perfectly understood. 

In a forthcoming paper, we will discuss more realistic f{R) gravity models with torsion confronting them with data. 
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